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We study a class of micromanipulation experiments, exemplified by the pulling apart of the two 
strands of double-stranded DNA (dsDNA). When the pulling force is increased to a critical value, 
an "unzipping" transition occurs. For random DNA sequences with short-ranged correlations, we 
obtain exact results for the number of monomers liberated and the specific heat, including the critical 
behavior at the transition. Related systems include a random heteropolymer pulled away from an 
adsorbing surface and a vortex line in a type 11 superconductor tilted away from a fragmented 
columnar defect. 
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Recent years have seen an explosion in the use of single 
molecule techniques to probe biological and other "soft" 
materials. It now possible, for example, to monitor the 
breaking of individual "lock and key" bonds [|| and the 
unfolding of individual proteins 10] ; the mechanical prop- 
erties of single DNA molecules [||| and the behavior of sin- 
gle molecular motors [Q have been characterized in great 
detail. In contrast to more traditional experiments, these 
new approaches give access to fluctuations on the scale of 
individual molecules, without the requirement for averag- 
ing over a macroscopic sample. One can, moreover, now 
push or pull directly on a micron-sized object, and watch 
how it responds. The potentially profound implications 
both for complex fluids and for biological physics — where 
single molecule techniques can often more closely mimic 
conditions in the cell than conventional assays — are only 
beginning to be explored. 

Despite a number of notable contributions, theory has 
often been out-paced by these rapid experimental ad- 
vances. Certainly, the tools available to analyze single- 
molecule experiments have not yet reached the level of 
sophistication and generality of theories of mesoscopic 
quantum systems. This is especially true when it comes 
to the role of quenched randomness, which, though often 
present, is typically neglected in initial theories of a given 
system. This Letter seeks to fill some of this gap. We 
study a class of micromanipulation experiments in which 
a polymer or other line-like object is pulled away from 
a confining potential well. An example of such a situ- 
ation is the pulling apart of the two strands of dsDNA 
(fig. |l|). Formally, the distance between the two strands 
may be viewed as the coordinate of a single polymer, 
and the base-pairing interactions between complemen- 
tary strands as a potential well. At a critical value of the 
pulling force, a novel phase transition occurs in which 
the two strands are pulled completely apart. Aspects of 
this transition for a homopolymer (or, equivalently, DNA 
with all base pairs the same) have been studied in a re- 
lated model of a flux line in a type II superconductor |^,^ . 
Here, we show that the transition is markedly different 
for random heteropolymers. In particular, the number of 
monomers liberated as the transition is approached di- 



verges much more strongly for heteropolymers than for 
homopolymers; similar differences appear in the specific 
heat. We calculate exact critical exponents and crossover 
functions for the random case. 

Figure |l| sketches the DNA-opening experiment: One 
of the two single strands of a dsDNA molecule is attached 
to a glass slide, while a constant force F directed away 
from the slide acts on the end of the other strand. Meth- 
ods for exerting a constant force on the piconewton scale 
have been developed by several groups Under the 

influence of the force the DNA partially "unzips" at 
one end, breaking m base pairs. In thermal equilibrium, 
the degree of opening m is of course a fluctuating quan- 
tity. Because the base sequence of protein-coding DNA 
appears to many statistical tests to be random and un- 
correlated along the backbone (at least up to a length 



M til 



Random 
DNA 

Sequence 



TO 




FIG. 1. Sketch of the "unzipping DNA" experiment. One 
of the single strands of a dsDNA molecule with a random 
base sequence is attached to a glass slide and the other is 
pulled away from the slide with a constant force F. As a 
result, the two strands partially separate, breaking m bonds 
(m = 2 in the figure). Inset: Schematic phase diagram in the 
temperature-pulling force (T-F) plane for a dsDNA molecule 
in 3 dimensions. At low T and F, the polymer is in the native, 
double-stranded state. At the phase transition line Fc{T), the 
DNA denatures, and the two strands separate. As indicated 
by the arrow, here we consider the unzipping transition that 
occurs when the transition is approached away from the F — 
melting temperature Tm. 
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scale set by the sequence's mosaic structure) the free 
energy landscape in which m fluctuates can be taken 
to have a quenched random component. Bockelmann, 
Essevaz-Roulet, and Heslot have performed an elegant 
series of experiments in a different statistical ensemble, 
measuring the average force required to hold the positions 
of both single strands fixed [^,0. However, because of 
subtleties associated with the thermodynamic limit in a 
single molecule system (see below), the two ensembles 
are not equivalent. 

In the remainder of this paper, we first introduce a 
coarse-grained model for the interaction of the two sin- 
gle strands of the dsDNA. By focusing on the unzipping 
transition induced by pulling on the single strands, we 
can avoid treating most of the degrees of freedom explic- 
itly, obtaining a problem that can be solved exactly by a 
mapping to a Markov process. Although for concreteness 
we will focus primarily on the DNA-opening realization 
of our model, our results also apply to a number of re- 
lated physical systems, some of which will be described 
in the conclusion. Throughout, we set fee = 1. 

The bulk melting transition of dsDNA (see inset 
to Fig. J^) can be described by a Peyrard-Bishop-like 
model |ll| . One views the two single strands as Gaussian 
polymers whose n^^ monomers have positions ri(n) and 
r2{n). Below the melting transition, it should be possi- 
ble to neglect non-native base pairings. The interactions 
between the two strands, coarse-grained over a number 
of bases, can then be described by a phenomenological 
potential energy term y„[r(n)] = [1 -\- fi{n)]h[r{n)\. Here 
his a, short-ranged attractive potential whose strength is 
temperature-dependent, r = ri — r2, and the variation 
with base sequence of the strength of the attraction be- 
tween strands is described by f}{n), which we take to be a 
random variable with short-ranged correlations. The ef- 
fective Hamiltonian then becomes, up to an uninteresting 
center of mass term, 

H.e. = f dn|™(£)%KHn)]} (1) 

where d is the spatial dimension and h is \/2 times 
the Kuhn length of single-stranded DNA. Standard 
arguments show that the partition function Zmoit = 
J 2?[r(n)] exp(— 7i,ncit/2^) obeys an (imaginary) time- 
dependent Schrodinger equation. 

The unzipping transition may be studied by adding to 
Hmcit the term T^puii = F ■ r(0) = F ■ r{N) - dn F ■ 
dr/dn. If the term proportional to r(N), which should 
not affect the opening near n ^ for a sufficiently long 
polymer, is dropped, a time-dependent version of the 
"non-Hermitian quantum mechanics" studied in |^ re- 
sults. In the time-independent case (corresponding to 
pulling apart homopolymeric dsDNA) there is a sharp 
first order unzipping transition at a critical value of the 
puUing force satisfying eo{T) = -F^^lP / {2dT), where 
£0 < is the ground state energy of the Hermitian quan- 
tum mechanics problem obtained by setting F = Q. In 



general, —F'^b'^/{2dT) is the free energy per monomer of 
the unzipped monomers aligned with the pulling force. 
The free energy per monomer of the dsDNA that has 
remained zipped is eo, independent of F. The physical 
interpretation of the unzipping transition is thus clear: 
For F < Fc{T), the DNA minimizes its free energy by 
remaining in the double-stranded form, while for F > Fq 
it is advantageous to pull apart as many bases as possi- 
ble. As F ^ i^c , the free energy difference between the 
bound, double-stranded phase and the pulled out, single- 
stranded phase becomes very small, and thermal fluctu- 
ations unbind a large number of monomers near the end 
of the DNA. As the transition is approached, the equi- 
librium ensemble average of the number m of monomers 
that are pulled out diverges like 

(m) - (Fc - F)"^ . (homopolymer) (2) 

Similarly, {{m— {m))'^) ^ (Fc — F)^^ near the transition. 
The thermal fluctuations about {m) are thus comparable 
to (to) itself. The divergence in (0) is analogous to the 
divergence in interface height near a wetting transition. 

We now determine how results such as (||) are mod- 
ified for a random DNA sequence. Sequence random- 
ness is at worst a marginal perturbation at the (F = 0) 
melting transition in 3 dimensions ||ri|-[l3||. The applica- 
tion of a Harris-like criterion however, shows that 
the same cannot be true for the unzipping transition: 
The typical variation per monomer in the base-pairing 
energy of a pulled out section of length < to > scales 
like < TO >~i/2^ y/Fc — F, which vanishes more slowly 
as the transition is approached than the relevant energy 
diff'erence |eo| - F'^b'^/{2dT) - Fc - F. To determine the 
correct critical behavior, we focus on the free energy cost 
of pulling out a given monomer. Define £{m) to be the 
free energy of a dsDNA molecule, subject to an applied 
force F, of which exactly the first m monomers are un- 
zipped. The change in £ from pulling out one additional 
monomer should have the form 

^ = / + r?(TO), (3) 
dm 

which may be integrated twice to obtain the partition 
function Z — dmexp[—£{m)/T] . Here / is the av- 
erage free energy difference between an unzipped and a 
bound pair of complementary monomers. It vanishes like 
Fc — F near the transition, and reduces to the familiar 
|eo| — F^fe^/ {2dT) in the absence of sequence randomness. 
The additional term 77 (to) takes account of sequence- 
dependent deviations from the average; it reflects the 
bare sequence (described by f}{m)), dressed by thermal 
fluctuations. As long as fj{m) is a random variable with 
only short-ranged correlations, it is reasonable to expect 
that ri{m) should also be short-range-correlated, with a 
correlation length on the order of the typical size f of the 
regions of local melting of the dsDNA strand On 
long enough scales, we can then take 77 to be Gaussian 
white noise, with correlator rj(m)rj{m') — A(5(to — to'). 
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where the overbar indicates a "disorder average" over 
the possible reaHzations of the quenched random base 
sequence. The parameters / and A may be calculated 
from the F = partition function ^mcit, for example in 
a low temperature expansion. The model summarized by 
Eq. (^) can also be derived from a discrete, Ising-like de- 
scription of the dsDNA ||l6| , and it still holds both when 
the single strands that have been liberated are character- 
ized as freely-jointed or worm-like chains and when there 
are significant excluded volume interactions 

The study of the unzipping transition can thus be re- 
duced to that of a single coordinate m in the random po- 
tential £{m). One immediate consequence is that there is 
no large parameter that defines a thermodynamic limit, 
and thus no equivalence between the ensemble consid- 
ered here and the conjugate ensemble in which m is held 
fixed. Below the unzipping transition, / > 0, and £{m) 
diverges with probability unity as m — > oo. In the en- 
semble studied here, the unzipping fork is thus always 
confined to the vicinity of to = 0. In the absence of ran- 
domness, the probability of unzipping to monomers is 
(//T) exp(— to//T), and one recovers, e.g. (^. If there 
is sequence randomness, the typical random contribution 
to £{m) is of order VAto; the random part thus exceeds 
the average contribution /to, which is responsible for the 



confinement, for 



< 



A// . This length scale diverges 



faster than 1// as / ^ 0, suggesting that a typical value 
of to might show a 1//^ divergence instead of the non- 
random 1//, with a crossover at / « A/T. 

Because of the confinement near to = 0, the unzip- 
ping transition does not exhibit self-averaging (see be- 
low). Nonetheless, one can still calculate averaged quan- 
tities and distributions over the possible realizations of 
randomness. To do this, one wishes to study the disorder- 
averaged free energy —Tin Z. The partition function 
Z{m) = dm' ex-p[—£{m') /T] of a finite-sized system 
of TO (bound or liberated) monomers satisfies 



dZ 
dm 



^^-£{rn)/T Z(0)=0; 



(4) 



Z follows simply by taking the limit of an infinitely long 
polymer: Z = limm^oo Z{m). Together, (||) and (^ form 
a system of coupled Langevin equations, analogous, for 
example, to those describing the Brownian motion of a 
massive particle, with £ playing the role of momentum 
and Z that of position. The associated Fokker-Planck 
equation for the joint distribution P{£, Z;rn} of £ and Z 
at "time" to follows in the usual manner flSl : 
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(5) 



By Laplace transforming with respect to Z and to to, one 
can solve (|^) and obtain an exact expression for the par- 
tial distribution J d£P(£, Z;m^ 00) in terms of modi- 
fied Bessel functions of order 2/T/A; —TlnZ and other 
thermodynamic quantities then follow by integration. 



The device of treating the quenched randomness as 
a Langevin noise thus leads to a number of exact re- 
sults. In particular, the average degree of opening (m) = 
-TdWZ/df satisfies 
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(6) 



where r'(z) = dr{z)/dz. The small / behavior of Eq. 
(|) is given by 

(to) ~ — ^ {Fc — F)^"^ , (random heteropolymer) 

(7) 

confirming our expectations for a crossover from a 1// to 
a power law when / ~ A/T. Similarly, the singular 
part of the heat capacity associated with the unzipping 
transition, C - d'^]^/dT'^, crosses over from a l//^ to 
a 1//"^ divergence. One can also compute the disorder- 
averaged values of higher cumulants of to. For small /, 
(to2) - (to) 2 = Td^h^/dp Unlike in the non- 

random case, the squa re root of this quantity diverges 
more slowly than (to) ^ 1//^, indicating that thermal 
fluctuations in to for a given realization of the quenched 
randomness (and for thus a given heteropolymer) typi- 
cally become small compared to the mean value as the 
transition is approached. 

A real space renormalization group approach to the 
model of equations (||) and (Q) due to Le Doussal, Mon- 
thus, and Fisher |l9| gives further insight into the unzip- 
ping transition. This technique gives leading order re- 
sults in the limit / ^ 0, where the authors have argued 
that it should be exact. In this limit, it allows one to cal- 
culate the distribution Q{{m)) of thermal average values 
over different realizations of randomness. This takes the 
form of a scaling function of (to)/2/A: 



ttA 



dw I 



(8) 



This distribution yields the same asymptotic behavior of 
(to) near the unzipping transition as the Fokker-Planck 

approach. It also predicts that ((to) — (to))2 ^ (to) 
for a polymer with a given random sequence of base pairs 
can thus deviate significantly from the disorder average, 
and this system is not self-averaging. 

For the randomness-dominated critical properties re- 
ported here to be observable, the variance A in the base- 
pairing energy must be sufhciently large. Then at the 
crossover from non-random to random behavior, / will 
also be large, and the typical value of (to) ~ T'^/A will 
be small enough that finite size effects do not become 
an issue. In this respect, dsDNA appears to be a very 
good candidate system. Under physiological conditions. 
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the difference in free energy of binding between poly- 
mers with only A-T base pairs and those with only G-C 
base pairs is of order T, meaning that (m) is only a few 
monomers when the crossover from pure to random be- 
havior occurs. 

In sum, we have described a randomness-dominated 
unzipping transition of dsDNA, obtaining exact expres- 
sions for the critical behavior and for the crossover from 
random to non-random scaling. Most notably, we find 
that when the base sequence is random and has only 
short-ranged correlations, the average degree of opening 
(m) diverges like l/(-F'c — F)'^ as the pulling force F ap- 
proaches a critical value Fc, in marked contrast to the 
l/(Fc — F) divergence found when all of the base pairs 
are identical. It should be possible to arrive at analogous 
results for the case of DNA whose base sequence has long- 
ranged correlations (as may be the case for non-coding 
DNA 1^). If a typical variation about the average energy 
grow s like , then balancing this against mf suggests 
(m) ~ 1//^/'^^^'^^; the short-range-correlated case is re- 
covered when P = 1/2. The biological significance of our 
results remains to be determined: Processes such as DNA 
replication and recombination often involve unzipping of 
the dsDNA. Usually, however, this is accomplished by a 
molecular motor relying on an outside energy source, so 
non-equilibrium effects must be considered. More gener- 
ally, the dynamics of the unzipping transition is an open 
question. 

We have focussed on the case of unzipping DNA, but 
our results hold equally well for a number of more conven- 
tional condensed matter systems described by the Hamil- 
tonian '^nieit + 'Wpuii • The pulling of a Gaussian ran- 
dom /leieropolymer away from an adsorbing surface is a 
natural extension of recent work on homopolymers [20|] . 
Other examples include a heteropolymer under tension 
pinned to a bulk defect [ pT|] , a magnetic flux line in a 
type II superconductor confined to a fragmented colum- 
nar pin and subject to a transverse field |]^,|l2|, and a 
simplified model of the corner wetting transition in two 
dimensions |^^. Related models are likely to be rele- 
vant to the transverse surface magnetization and surface 
specific heat near the Bose glass transition of a bulk su- 
perconductor 1^ and to adhesion in a random environ- 
ment ||. 

After this work was submitted for publication, we 
learned that related results had been obtained, in a dif- 
ferent physical context, for a discrete version of Eqns. 
(H) and jH). It is a pleasure to thank D. Branton, 
D.S. Fisher, and T. Hwa for helpful conversations and T. 
Hwa for bringing to our attention. This research was 
supported by the NSF through grant DMR97-14725 and 
through the Harvard MRSEC via grant DMR98-09363. 
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